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VELOCITY DISTRIBUTIONS ON TWO-DIMENSIONAL WING-DUCT INLETS 

BY CONFORMAL MAPPING 

By W. Perl and H. E. Moses 


SUMMARY 

The conformal-mapping method of the Cartesian mapping 
function is applied to the determination of the velocity distri- 
bution on arbitary two-dimensional duct-inlet shapes such as 
are used in wing installations. An idealized form of the actual 
wing-duct inlet is analyzed. The effects of leading-edge stagger, 
inlet-velocity ratio, and section lift coefficient on the velocity dis- 
tribution are included in the analysis. Numerical examples 
are given and, in part, compared with experimental data. 

INTRODUCTION 

Inlet contours for wing-duct installations, such, as those 
used to conduct cooling air to engines, are generally designed 
on a more empirical basis than airfoil sections because the 
geometry of a wing-duct inlet, and hence the determination 
of its velocity distribution, is more complex than that of an 
airfoil section. By means of the eonformal-mapping'method 
of reference 1, however, the ideal incompressible velocity 
distribution over two-dimensional wing-duct inlets for arbi- 
trary lift coefficients can be calculated with about the same 
labor as in the corresponding calculation for an isolated 
airfoil. 

This method was applied to an arbitrary two-dimensional 
wing-duct-inlet section at the NACA Cleveland laboratory 
in 1945 and the application is presented herein. The theory 
is illustrated by numerical examples, which are, in part, 
compared with experimental data. 

ANALYSIS 

SYMBOLS 

The more important symbols used in the paper are listed 
here. All velocities are expressed as fractions of the free- 
stream velocity; that is, the free-stream velocity is taken 
as unity. 

c chord of duct-inlet section 

Ci section lift coefficient 

d horizontal distance between leading edges of duct 

inlet 

h vertical distance between leading edges of duct 

inlet 

r stagger ratio in f -plane 

s stagger ratio d/h in z-plane 


v velocity on surface of duct inlet 

v D<0 velocity infinitely far inside duct 
i’a velocity at duct inlet 

p plane of circle 

z plane of duct-inlet section (x-fiy) 

f plane of chord lines (£+hj) 




(b) -Analyzed. 

Figube 1.— Wmg-duct Installation. 


THE CONFORMAL TRANSFORMATION 

The actual two-dimensional wing-duct configuration of 
figure 1(a) is replaced by the contour shown in figure 1(b). 
The two changes made in the original configuration are (a) 
removal of the internal flow resistance, and (b) replacement 
of the streamlines A'B' and C'D' by the parallel, straight, 
rigid boundaries AB and CD. Change (a) results in a flow 
field of constant total pressure, and change (b) in a simply 
connected flow field. The analysis is thereby considerably 
simplified. Both effects associated with the replaced features, 
namely, variable inlet-velocity ratio and angle of attack, 
respectively, can be adequately represented in the flow 
function for the simplified configuration. For conventional 
wing-duct installations, the region of interest at the inlet, as 
regards velocity distribution, is sufficiently far from the 
region in which changes (a) and (b) were made that their 
influence on the required velocity distributions is negligible. 
(See section Illustrative Examples.) 
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The simplified duct-inlet contour in the 2 -plane is now con- 
formally mapped onto the staggered semi-infinite parallel 
straight lines, QAB and PCD, in the f-plane (fig. 1 (b)). 
This mapping is accomplished by the Cartesian mapping 
function (CMF), defined as the vector distance 2 — f between 
conformally corresponding points in the 2 - and f-planes 
(reference 1); thus 


2— r =(*—£) +i(y—v)] 
—Ax+iAy \ 


The calculation of the CMF is carried out by considering it 
as a function of the central angle 4> of the p-plane circle into 
which the f-plane contour can be conformally mapped by a 
known transformation. Inasmuch as z— f is regular on and 
outside of the 2 - or ("-plane contours, by the conformal 
transformation from f to p it is also regular on and outside 
of the p-plane circle. The real and imaginary parts of the 
CMF on the circle itself are therefore related by 


Ax (4>) = - -^j^Ay (<*>') cot ~~ d<t>' (2) 


Ay(4>) =-^j^Ax(<t>') cot d<t>' (3) 

The conformal transformation of the f-plane, staggered 
semi-infinite parallel lines, into the p-plane circle is carried 
out in two steps. In the first step a Scliwarz-Christoffel 
transformation takes the f-plane polygon into the real axis 
of a f-plane such that the upper-half f-plane corresponds to 
the f-plane. With the correspondence _of_boundary points 
indicated in figure 2, this transformation is (reference 2): 





(ws+tfOf+ttsMt log, fJ+Ci (5) 

The six constants given by u%, (real), and C u C z (complex) 
are determined for the orientation and the scale indicated in 
figure 2 by the six conditions: 

(a) Ci real (staggered lines horizontal) 

(b) a=l, scale factor in f-plane 

(c) r equals desired stagger in f-plane, f(uj)=r-|-r— i 

(d) «j=“- 1, scale factor in f-plane 

(e) upper leading edge in f-plane at point (r, 0) or 

f (w 4 ) = r (two conditions) 

The constant r is inserted in condition (e) in order to 
locate the leading edge of the upper inlet section tangent 
to the y-axis. By use of the foregoing conditions, equation (5) 
reduces to 

f-i (s- 0 [1 (s- >)+ log - s +T (6) 

The quantity m=Ut is the following function of the stagger 
ratio r: 

7rr=log, (™“) W 

Equation (7) is plotted in figure 3. 



Stagger in $-ptane, r 

F tor be 3.— Stagger constant as a function of stagger In {■-plane. 
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The second step of the desired transformation from f to p 
consists in mapping the upper-half f-plane onto the region 
outside of the p-plane unit circle "by a bilinear transformation, 
here taken as 


t=i 



( 8 ) 


The correspondence of points for equation (8) is indicated 
in figure 2. The use of other bilinear transformations is 
discussed in the section Illustrative Examples. 

Equations (1), (6), and (8) constitute the conformal trans- 
formation from the region around the duct-inlet section in 
the physical 2-plane to the region outside the unit circle 23 - 
plane. These equations, with give for conformally 

corresponding points on the boundaries 


f=cot 


± 

2 


(9) 


x=£+Ax 

*=Km COt i -1 )S(= COt 5 -1 ) +1 ] - 


t (10) 


1 1 1 


cot ^ +*A a*(<£)+r 


y=ij+Ay 

y=Ay{$) 

y=Ay(4>)-l 


0 OO upper-duct inlet section >- 
r<<£<2r lower-duct inlet section J 


( 11 ) 


The leading-edge points of the upper-duct- and Iower- 
duct-inlet sections may be defined as the upstream points of 
tangency of normals from the “chord” lines OA and RC 
with the duct-inlet contours (fig. 1 (b)). These points may 
be found as functions of <f> by minimizing x with respect to 4> 
in equation (10). The resulting condition is 


dAx (cot|-m)(cotf+l) 
d<j> itm sin $ 


VELOCITY DISTRIBUTION 


The velocity distribution on the duct-inlet section is given 


by 


v= 


dW\ 

dz 


in which the complex potential ll T is 

W=?+- f+? log, t 

X 7T 


(13) 


(14) 


The term f represents a uniform flow velocity to the right, of 
unit magnitude in the f-plane, and gives a free-stream veloc- 

ity of unity in the physical plane. The term — t represents 

a uniform flow in the t-plane and corresponds to a circulatory 
flow around the duct inlet in the physical plane. This 
term gives the effect of angle of attack on the physical duct- 
inlet section, although the geometric angle of attack of the 
section analyzed must remain zero because of its semi- 

jB 

infinite extent. The term — log, t represents the flow due 

T 


to a source at the origin in the f-plane and gives the desired 
inlet velocity into the duct in the physical plane. 

The quantitative effect of the parameters A and B in the 
physical plane is determined by evaluation of the complex 

velocity 

dlI T f ,A dt.B dt\ df . . 

n = v + ^d^vtWTz - (15) 

where, by equation (6), 


dt mnt 
df~(<— m)(f+l) 


(16) 


and, because z— f is regular on and outside of the p-plane 
circle, 




07) 


dz_ dp ^ 7ic„ 
.2f -1 


dz . 2i xmf 

df _1 + (f-i) 2 (f-m)(f+l) 



(18) 


Infinitely far inside the duct in the physical plane, the corre- 
spondence of points is: z= c°, f = m by equation (17), t—0 
by equation (6), and p=—l by equation (8). . Hence, at this 

point, ^=1 by equation (18), ^=0, and y ^=— irby equa- 
tion (16); and equation (15) gives for the velocity v Dm infi- 
nitely far inside the duct 

v Da ^l-B (19) 


The velocity distribution on the inner wall of the duct-inlet 
section becomes almost constant a short distance behind the 
leading edge. (See section Illustrative Examples.) The inlet 
velocity v u is defined as this asymptotic value. The inlet ve- 
locity t’ a will be different from v Da if the height at the inlet is 
different from the height (unity) Infinitely far inside the duct. 
Infinitely far upstream of the duct-inlet section the correspond- 
ence of points is: z= — f=— °°, f=i°°, and p= 1, and 

l,^j;=i^=0. This result holds infinitely 

far outside the duct in any direction. Hence, the free-stream 
velocity is by equation (15), unity. 

The quantity A may be evaluated either as a function of 
the stagnation-point locations on the duct inlet or as a function 
of a suitably defined lift coefficient. In terms of the stagna- 
dW 

tion points, given by -g— • =0 in equation (15), and with equa- 


consequently 


tions (16) and (9), 


(m— 1) cot 
A= 


‘I'-cot* 
m cot 


,( 20 ) 



-[l + m(^-l)]±V[l+w(H-l)] 2 +4m(l-)5) 
2 


( 21 ) 


For a given A (and B ) equation (21) is a quadratic equation^ 
for the two stagnation-point locations. When quantity A is 
alternatively regarded as a function of lift coefficient e t , the 
section lift coefficient is defined in terms of circulation and 
chord by the well-known isolated-airfoil relation 
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C=j - - - ( 22 ) 

The cliord c is defined as the over-all. length of the. wing- 
duct-inlet section in the free-stream direction (OE in 
fig. 1(b)), and the circulation r, as the line integral of the 
velocity over the circuit CFGHKMAEC around the wing- 
duct installation. This circulation can be evaluated as the 
sum of the potential difference over. the. lower surface $ G — $ c , 
and the potential difference over the upper surface $ E — $ H . 
The difference of potential over the paths GH and EC is 
neglected because the velocity is here approximately per- 
pendicular to the path. Hence, by equation (14) _ 

$ E — $h=?e — (^ a £h)H log* ~ 

ir tt j H 


$G — $C = £g~€c+“~ (fa - •fcl+— loge ~ 

IT IT £ c 

and 


r=(4 >E“ (^g — $ c) 


r= (£e~ Zc) — (£h — £g)'+~ — fc) — (fn— f g)J + L 



(23) 


Finally, when equation (23) is solved for A, and T is expressed 
in terms of Cj by equation (22) with c=z E , 


PROCEDURE FOR CALCULATION OF CMF 

The calculation of the CMF A x(<j>), Ay(4>), and stagger 
constant m for a given duct-inlet section may be carried out 
by a process of successive approximations similar to that of 
reference 1. The steps are outlined as follows: 

1. The duct-inlet section is drawn in normal form 
(fig. 1 (b)). Point 0 is the origin and the scale is such that 
the normal distance between the chord lines OA and RC is 
unity. The stagger s=djh of the duct inlet is, in general, 
different from the stagger r of the chord lines. 

2. A set of abscissas x(4>) is calculated for a standard set 
of values of <j> by equation (10). The A x(<j>) and r may be 
that of_a previous example or, at worst, equal to zero. .The 
value qI.m may be taken from figure 3 for r=s. 

3. The. ordinates y of the duct-inlet contour corresponding 
to the abscissas x of step 2 are measured. The function 
A y(<£) is thereby determined (equation (11)). 

4. The function A z(<j>) is calculated from Ay(<f>) by equa- 
tion (2). 

5. The. functions Ax(<f>) of step 4, A y{<f>) of step 3, and m 
of step 2 constitute by equations (10) and (11) a duct-inlet 
section' of which the difference in abscissas between the 
leading edges is, in general, other than that specified. The 
constant m is therefore adjusted to make this difference 
equal to the specified value. To this end equation (12) 
(corresponding to the values 0, and 0 2 for the two extremities) 
and the equation for the difference d in the leading-edge 
abscissas. 

x(m, <£i) — x(m, (fe) —d (27) 


it ^sse— (Ie— £c) + (£h— £g) — B log, (^0 


( < e — tc )— (t h —£g ) 


(24) 


The quantities x, £, and t at the various points indicated in 
equation (24) are given in terms of the corresponding central 
angle 0 by equations (9) and (10). The various. 0 values 
are known when the conformal transformation of the duct 
inlet has been carried out. 

For a CMF A r(0), A y(4>), stagger constant rrij and the 
constants A and B corresponding to the lift coefficient and 
the inlet-velocity ratio, the velocity distribution on the 
duct-inlet contour is given by the absolute magnitude of 
dWjdz on the boundary. On the boundary, p=e*+, 


dz ., , d{z— £) ^ 

3r 1+ “ar" '• 

dz_ d{Ax+iAy) d<j> dt 
dl; 30 3i 3| > 


(25) 


obtained from equation (10), can be solved simultaneously 
for 0i, 02, and m. A more convenient procedure is one of 
iteration. Initial values 0! and <j> 2 for minimum x arc 
graphically obtained by plotting equation (10) in the neces- 
sary regions. A value of mis then obtained from equation (27). 
With this value of m, values of 0i and 0* are again 
graphically found for minimum x by equation (10). The 
process.is continued until 0 lf 0 2 , and m do not .change ap- 
preciably in successive calculations. Finally, a constant r 
is so chosen that z(0i)=O. The derived inlet section is 
now in normal form. 

6. The values of m and t derived in step 5 and A x(4>) and 
A y(0) of steps 4 and 3 yield a shape by equations (10) and 
(11), which can be compared with the given one. If the 
agreement is not sufficiently close, steps 3 to 5 are repealed. 

7. The velocity distribution is obtained by substitution 


of the final m and the derivatives 


dAx 

14 




of the final 


Substitution of equations (9), (16), and (25) inequation (15) 
yields for the velocity distribution on the duct-inlet section 


(cot t—m) (cot ^+l\+mA cot 


J 

(cot | -m) (cot |+ 1^) 

dAx 

2 

+1 

( dAy'A 

'V 

Ttm sin 0 

d<f> _ 

\d4> ) 


( 26 ) 


CMF in equation (26). The value of B is chosen to produco 
the desired inlet velocity (the velocity given by equation (26) 
on the inside walls of the duct-inlet section). The value 
of A. Js_chosen to locale the stagnation points in the 
desired manner (equation (21)) or for a desired nominal lift 
coefficient (equation (24)). 

The inverse problem, namely, the calculation of the duct- 
inlet section to produce a prescribed velocity distribution, 
may be treated by the methods given in references 1 and 3. 
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ILLUSTRATIVE EXAMPLES 

As a first application of the theory, the symmetrical wing- 
duct installation (m=1.0), on which pressure distributions 
were measured in reference 4 (shape 9), was analyzed. The 
installation is shown in figure 4 and the ordinates are listed 
in table I. The trailing-edge portions were actually flaps 



Figure 4.— Symmetrical duct -inlet section (reference 4, shape 6). m-1.0. 


by which the inlet velocity was varied. The scale used for 
the calculation was such that the distance between trailing 
edges was unity, as assumed in the theory. An evenly 
spaced set of 48 ^-values was taken of which only 24 were 
actually used because of the symmetry. Of these 24 values, 
21 were included in the front 8 percent of the chord. This 
portion of the duct inlet was therefore the portion effectively 
analyzed. The leading-edge portion is plotted in figure 5 
to a scale such that the vertical distance between the leading 
edges, the entrance height h, is unity. 



Figure 5.— Leading edge of symmetrical duct-inlet section with t=0 and m— 1J3. 


The CMF obtained after four approximations (which 
produced coincidence of the specified shape and the derived 
shape) is listed in table II and plotted in figure 6. In the 



first two approximations, the airfoil was drawn to an abscissa 
scale of 25 inches for the chord and had an ordinate scale 
four times the abscissa scale. The last two approximations 
were made for the airfoil drawn to a scale such that the 
chord length was 100 inches; the ordinate scale was the same 
as the abscissa scale. The values of Ax were computed, for 
the most part, by the method of numerical evaluation of 
conjugate functions developed in appendix C of reference 5. 
Near 0° and 180°, because of the rapid variation of A y(<f>) 
in these regions (fig. 6), Ax was obtained by plotting the 
integrand of equation (2) and graphical integration. The 
values of the CMF graphically obtained are indicated in 
table II. The velocity distributions, also listed in table II, 
were calculated for inlet-velocity ratios v K of 0, 0.5, and 1.0 
and for nominal lift coefficients of 0, 0.3, 0.6, 0.9, and are 
shown in figure 7. The derivatives of the CMF used in 
calculating the velocity distribution were obtained by graph- 
ical measurement from the CMF. 

The velocity distribution for cj=0 and t*»=0.5 .satis- 
factorily checked that experimentally obtained in reference 4 
for c { =0 and i? a =0.473 (fig. 7(c)). The reason for the 
discrepancy between theoretical and experimental inlet- 
velocity ratios at which the velocity distributions. agreed is 
not clear. Possible reasons are the changes in downstream 
shape required by the analysis and a difference in the method 
of specification of inlet-velocity ratio. The theoretical inlet 
velocity v n has been defined as the constant value approached 



Velocity, v Velocity, 
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Velocity, u __ _ ^ ^ ^ Velocity, 
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Figure 7.— Concluded. 


by the velocity on the wall of the duct at a short distance 
behind the leading edge. The experimental determination 
of the inlet-velocity ratio in reference 4 was not made 
entirely clear. 

The velocity distributions for the c t values 0.3 and 0.6 
were also compared with the experimental data of reference 6 
obtained for the same duct-inlet section at various angles of 
attack at a Mach number of 0.20. The comparison (given 
in figs. 7 (c) and 7 (d)), indicates the validity of the theoretical 
analysis, particularly of the derivation of the nominal section 
lift coefficient c t . 

The feature of the velocity distribution shown in figure 7 
that should be particularly noted is the closeness to the 
leading edge (well within the 8-percent of the chord length 
that was studied) at which the greatest changes in velocity 
distribution occur as a result of a change in operating con- 
ditions r* or c s . This fact justifies and requires the analysis 
of a region very close to the inlet, that is, the concentration 
of the chosen set of ^-points close to the inlet. 

In order to illustrate the use of the theory for the staggered 
case, m^l.O, the CMF Ar(0) and A y(4>) for the symmetrical 
inlet was used with the m-values 1.5 and 2.0. These shapes 
and velocities are shown in figures 8 to 11 and are given in 
tables III and IV, respectively. In the graplis of the duct 
inlets, the ordinates have been so adjusted that the upper and 
lower leading edges are at 0.5 and —0.5, respectively. Al- 
though the derived shapes are different from the original 
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Fioubk 9. — Velocity distribution on upper and lower leading edges of ndnsymmetrlcal 
duct-inlet section with *-0.132 and y»-0.5. 


unstaggercd one, evidently tlie effect of stagger is to reduce 
the velocity peaks for positive lift coefficients. 

When a more highly staggered inlet was derived for 
m=3 r .O by the foregoing method, the upper contour of the 
resulting inlet was found to be excessively thick. The points in 
the physical plane corresponding to A x(,<f>), Ai/($), and m=3.0 
were therefore rearranged by using the same A/ and Ay, 
regarded, however, as functions of 6, with 9 related to <i> 
by the bilinear transformation (see appendix B of reference 5) 


P 


pH 

n— 1 
it+1 


n— 1 
71+1 

PHI 


(28) 


in which, on the p- and p'-pl ane unit circles corresponding 
to the duct-inlet contour, 

—e tf (29) 

The choice n— 1.5 produced the shape shown in figure 12. 
The ordinates are listed in table V. It should be noted that 
the use of an auxiliary bilinear transformation (equation (28), 
for example) provides a very flexible and convenient 
method of distributing a given number of mapping points 
in the optimum manner. The auxiliary bilinear trans- 
formation may also be used to smooth out a sharply peaked 
function of 4> to make its conjugate more easily calculable. 
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Fiqcee 11.- —Velocity distribution on upper and lower leading edges of nonsymmetrleal 
duct-inlet section with t— 0.215 and 77»— 2J3. p«=0.5. 



Fiqcee 12.— Leading edge of nonsymmetrleal duct -Inlet section with #=0.280 and m=3.0. 

•-US, 

i 

CONFORMAL MAPPING OF LEADING-EDGE REGIONS 

The requirement that the velocity distribution need be 
accurately known only near the leading edges permitted the 
great simplification in the mapping consisting ip. replacement 
of the doubly connected region by a simply connected region. 
The modification of the contour shape far behind the leading 
edge did not appreciably alter the velocity distribution at the 
leading edges. Corresponding simplifications can be effected 
in other problems involving conformal mapping of aerody- 
namic shapes where only the leading-edge region is of interest. 

Thus, for example, the leading-edge region of an isolated 
airfoil can be regarded as joining a semi -infini te shape, as 
indicated in figure 13. The mapping of such a contour into 
a circle is quite simple. The leading-edge contour, 2 -plane, 
is mapped onto a semi-infinite chord line, f-plane, by the 
CMF 

z—£=Ax+iAy= (x—g) +i(y—v) (30) 

The semi-infinite chord line is mapped onto an infinite 
straight line, t-plane, by 

f=£ 2 (31) 




78 


HE PORT NO. 893 — NATIONAL ADVISORY COMMITTEE FOR 'AERONAUTICS 



Figure 13.— Mapping of leading-edge region fbr isolated airfoils. 


and, in turn, the t-plane contour is mapped onto a unit circle 
by a bilinear transformation, such as 

' • (?'!) < 32 > 

On the unit circle p—e { *, equations (30) to (32) give for the 
coordinates x and y of theleading-edge contour in the physical 
plane 

x=cot*| +Ax(<j>) (33) 


SAW-TOOTH FUNCTION AS INITIAL APPROXIMATION 

In "the mapping of semi-infinite contours, it may be re- 
quired, or it may be simpler, to consider a contour for which 
the thickness at infinity is finite, as indicated in figure 13. 
The ordinate function A y(<f>) will in this case be discontinuous 
at the value of <j> corresponding to the point at infinity on the 
physical-plane contour. The calculation of the CMF for 
such a contour will be simplified if the Cartesian mapping 
function Ar.+iAy is considered as the sum of two component 
Cartesian mapping functions AiX+iA^ and AjX+tAjj/, of 
which Atf+iAiy is analytically known and represents a 
contour with the same thickness at infinity. 

Thus, if the thickness at infinity of the contour in the 
physical plane is T, the “first harmonic” saw-tooth function 

(%• M 

AtfM-l'fJ-f) (37) 

will yield a shape with this thickness. The function Aix(<j>) 
conjugate to A } y(<j>) maybe simply obtained from the integral 
relation for the conj'ugate derivative (equation (C3) of 
reference 5) 


dAz(c£)_ 1_ p*- 

d<t> 4ir Jo 



d<f>' 


(38) 


Substitution of equation (37) in equation (38) and integration 
(using integration by parts) yields 


dA x x T . <(> 

d<f> 2 ir cot ’ 2 


y=Ay(<j>) 


(34) 


which by integration gives for AiX 


The mapping of the leading-edge contour by equations (33) 
and (34) involves little more than the calculation of 
conjugates. 

The velocity distribution is obtained from the complex 
potential 

W=J+At (35) 

in which the term £ represents the uniform free-stream flow 
and the term At a circulatory flow around the leading edge. 
On the leading-edge contour . the velocity distribution 
| dW/dz{ becomes . 



A similar development can be made for a cascade of leading- 
edge regions, which bears the same relation to a cascade 
of airfoils as the leading-edge region just treated bears to 
the isolated airfoil. 


T <t> 

Aix= — log, sin | (40) 


The ordinate function derivative is evidently 



Figure 14.— Saw-tooth ordinate function*. 
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The CMF for a “second harmonic’' saw tooth (fig. 14) 
corresponding- to a duct-inlet section with thickness T (of 
each component contour) at infinity may be obtained from 
that of the first harmonic saw-tooth function by replacing 
<j>l 2 in equations (37), (39), and (40) by <j> and by doubling 
the derivatives. Thus, for the second harmonic saw-tooth 
function 


&iV = ~ (f —$) ooo 

r 

Ai3/=~(^y — $) tO<2s-^ 


(42) 


d&iy = T 

d<j> x 


T 

A. l x=— log, sin <f> 

* TC 


(44) 


dAi x_T 

d<f> x 


cot <f> 


(45) 


The duct-inlet shapes corresponding to the second harmonic 
saw-tooth CMF have been calculated by equations (10) and 
(11) with m=1.0 (no stagger) and for T=0.1, 0.2, and 0.3. 
The contours are shown in figure 15. For these shapes, the 
velocity infinitely far inside the duct is, by equations (43), 
(45), and (26) with <0=r, 



F lgf ec 15.— Leading-edge shapes for duct inlet using the second harmonic saw-tooth function. 

a-0; ot-L 


Vd co — 


1-5 

i-t 


(46) 


which, when compared with equation (19), shows the effect 
of narrowing the duct at infinity. 


Aircraft Engine Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 1, 1947. 
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TABLE I— ORDINATES OF SYMMETRICAL DUCT INLET 


[Reference 4, Bhape 8] 



Ordinate of 

Ordinate 

Station 

(percent 

chord 

from 

leading 

edge) 

outer 

surface 

(measured 

offnner 

surface 

(measured 

from center 

from center 

line of 
channel) 
(percent 

line of 
channel) 
(percent 


chord) 

chord) 

0 

8.343 

3.34 

.25 

8.060 

3.07 

.5 

3.835 

3.05 

.75 

3.876 

3,05 

L 25 

4.228 

3.07 

2.5 

4.745 

3.13 

5 

5.532 

3.25 

7.5 

6.137 

3.36 

10 

0.652 

3.50 

15 

7.467 

8.78 

20 

8.088 

4.08 

! 25 

8.693 


30 

8.866 

4.67 

85 

8.224 


40 

9.370 

4.90 

45 

50 

0.435 

8.391 

6.10 

55 

8.240 


60 

8.666 

4.94 

65 

8.610 


70 

T 1 

4.27 

76 





2.54 






L55 





1.314 

1.31 






00 

o 


TABLE II— ORDINATES AND VELOCITY DISTRIBUTIONS FOR SYMMETRICAL DUCT INLET WITH s=0 AND m= 1.0 

I [r-0; **,-76.00°; <6 e -286.00“; * K -0.882°; * c =363.12°; A-24800] 


<A 

(deg) 

x [h 

m 

H 

AX 

Ay 

Velocity, c 

(o.-Oj B-1.00000) 

(n.-OJS; £--0.38907) 


ei-0 

0.3 

0.08494 

0.0 

1.16989 

0.9 

1.75483 

0 

0 

0.3 

0-58494 

0.6 

1.16888 

0.9 

1,75483 

0 

0 

03 

0,58494 

0.0 

1.16989 

0.9 

( 1.75483 

Upper duct Inlet section 



0X7. S': 

2 

3 

4 
6 
6 

7 

8 
9- 

10 
11 , 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 

. 22 

23 

24 

12,3640 
'2. 8117 
1.0887 
.6162 
.2490 
' ; 1226 
,0624 
■„ ''. 0221 
. '.0041 
.0017 
, . 0038 
' ■ . 0130 
•: i .0269 
’ ■ i .0471 
.0718 
' ;:.1034 
1 -.1386 
.1842 
.2363 
.3067 
.. 1.3912 
1.6227 
.7401 

CO 

ifcJi. . 

“I*'' 

1, 2098 
,9259 
.7713 
.6696 
6121 
.5708 
.5410 
.5166 
.4996 
.4858 
.4766 
. .4713 . 
.4690 
.4682 
.4676 
.4690 
.4705 
.4786 
,4766 : 
J4789 
.4843 
.4950 ' 
.2011 
* •> ,4 

i •—2, 8592 
•—6. 1796 
— 1.2696 
—.5144 
-.2265 
1300 
-.0546 
-.0115 
.0400 
.0099 
.1000 
.1307 
.1623 
,1879 
.2268 
1 .2655 

.3168 
.3685 
.4379 
.5131 
. .6195 

.7400 
.9389 
•1.2595 
• 3,0449 

’1.0534 
- 2.6076 
" 1.8018 
1.4174 
1.1644 
L 0217 
.9190 
• : 78448 
.7839 
.7420 
.7078 
' .6849 
, .67?0 
.6659 
.6640 
.6621 
..6659 
.6697 
. .6773 

.6849 
•' .6906 
. ' ;7040 
:7206.. 

0 

l, 0000 V 
1.0697 
1.1719 
1. 1944 
L 2238 
L 2535 
L 8460, 
L4799 
1.6354 : 
1.8827 
2.0239 
1.7745 
1.2056 ' 
.7651 
,.4713 : 
.2991 
.1961 
. 1305 ‘ 
,0858 
.0543 
.0324 
a 0172 
.0069 
.0017 

1.U07 
’ L 2622 
1.3334 
L 4157 
1.5025 
0721 
1.9068 
2: 1877 
2.6185 
2.9324 
2.6848 
1. 9109 
1,2754 
.'8305 
.5609 
.3947 
.2853 
.2069 
» 1480 
. .0319 
. .0677 
.0374 
.0167 

0 

1.0000 
1. 1617 
1.8624 
1.4724 
1.6076 
1. 7514 
1.9982 
2.3337 
2.7400 
3.3544 
3.8408 
8.5951 
2.6101 
1.7857 
1. 1898 
.8227 
.5938 
.4401 
.3280 
.2416 
.1740 
.1182 
.0079 
.0318 

0 

.1.0000 . 
1.1927 
1.4427 
1.6113 
1.7993 
2.0003 
2.3243 
2.7606 
3.2923 
40902 
47493 
45063 
8.3213 
2.296Q 
1.6490 
1.0845 
.7920 
.5948 
.4491 
.3353 
.2447 
.1687 . 
.0984 
.0468 

0 

n.0000- 

'1. 0638 
-^.1437 
him 
1,1018 
1.0529 
ti0252 
■ .9800 
. .8782 1 
.7161 
.3686 
.1212 
.4691 
, |V 0167 
■.6406 
.6314 
.6210 
.6148 
.6086 
.6008 
.5948 
.6859 
, 5435 
.5432 
ly. 3891 

f 1. 2340 • 
1,2678 
L 2936 

Zzow 

- 1.3513 , 
1.4069 
V 1.4305 
% 1,4509 

I. 2770 

• .7891 
-2301 
* 1064 

.. - 28X3) 
';369fl 
.4223 
. ,4601 

• , ,4875 

.8071 
, .5240 
•. ,5854, 

■j .6130, 

. 5281 1 

J, 389i; 

j : i 'i 

t'3242 
,t:4O07 
1*4854 • 
1/5508 
U«774 
1. 8338 • 
1.9828 
2.1868 
2.1855 
1. 0993 
. 9414 , 

. . 4039 , 

;0780 
id078- - 
>2237 1 
13053 
.'3664 
.4135 
.4532 
.4849. 
.4825 
.6181, 
1.3891': 

i 1 >■' 

•• i>0000 
11863 
1/4144 
L.6467 
l* 0773 
L7997 
'2.0035 
2,2607 

- 2.6361 
2.9226 
3,0940 
2.6096 
. 1. 6467 
, ,9142 
' v4372 

.1540 
.*■ 0250 
■1.1505 
.2463 
.3198 
v . 3825 
.... .4844 
J4519 
.4980 
1.3891 

- •, :■ 

1.0000; 
1.0680 
1. 1201 
1.0739 
.9998: 
. .8862 
.7671 
.6022 
.2154 
.2807 
1. 0160 
1 7066 
1,8087 
1.7710, 
1.8597 
1.4090 
1.3038 
1.2378 
1.1889 
1.1483 
1.1190 
1090Q 
1.0034 
.9986 

Ltogq 
1. dwe- 
lt 2104 
1.2129 
1. 1916 
1, 1841 
1. 0832 
.9891 
i7977 
.4751 ' 
.1065 
.7962 
1635 
1.2013 
1.2105 
1.1472 • 
1.1051 
1.0836 
L 0678 
L 0546 
1.0482 
1.0395 
.9729 
.9835 

; fc?400 
L 3006 
1.3518 
1.3834 
1.8830 
1.4093 
' 1.41Q0 
1.3500 
1.2110 
.8020 
,1151 
,4583 
,7510 
.8512 
,8854 
.0065 
.9282 
.9467 
.9610 
.9774 
. >9890 
*9424 
.9684 

i.l.OQOO 
1. 1810 
1. 3009 
■1. 4608 
L 0752 
1. 0319 
X. 7364 
1.8429 
■ 1.9023 
' 1.9408 
1.7104 
1.0263 
.2470 
.2407 
. .4920 
.0230 
.7079 
.. ,7734 
.8200 
,8673 
- .9007 
.9386 
: .9119 
.9684 

lower duct-inlet section 


1 

o. 7m 

.5227 
-8913 
.8067 
.2363 
.1842 
.1885 
.1034 
,0718 
,0471 
.0259 
.0136 
.0033 
.0017 
.0041 
.0221 
.0524 
.1225 
.2490 
.5162 
L 0687 
2.8117 
13.8540 

-a 2011 

—.4960 
— . 4843 
-i 4789 
-,4706 
-.4736 
-,4706 
—i 4G90 
-.4676 
-.4682 
-.4690 
-. 4713 
-.4766 
-.4868 
4996 
5165 
-,5410 
-.5708 
6121 
-.6695 
-. 7713 
-.9259 
-1.2096 
-.8662 

• 8.0449 

• 1. 2595 

.9389 
.7400 
.0195 
.5181, 
,4379 
.3685 
.8168 
.2655 
.2263 
.1879 
.1023 
.1307 
.1060 
.0699 
.0400 
-.0115 
-.0546 
-.1306 
-.2365 
-.5144 
— L 2595 
•-S 1796 

0 

-.7306 
-.7040 
-.6900 
-. 6849 
-.6773 
—.6697 
-.6659 
-.6621 
-.0640 
-.6659 
-.6716 
-.6849 
-.7078 
-.7420 
-.7839 
-.8448 
-.9190 
-1,0217 
-1,1644 
-1.4174 
—1.8018 
-2,6070 
— leOfiSt 

0 

.0017 
.0069 
.0172 
.0324 
,0543 
.0858 
.1305 
.1961 
.2991 
.4713 
,7651 
1.2066 
L 7745 
2.0239 
1.8827 
1.6354 
L 4799 
1.3460 
1,2535 
1,2238 
L 1944 
b 1719 
1.0607 

0 

.0134 
.0237 
.0333 
.0384 
,0393 
.0854 
.0248 
.0026 
.0375 
.1120 
.2548 
.5004 
.8642 
1.1155 
1.1465 
1-0881 
1.0631 
L0198 
L 0046 
1.0320 
1.0554 
L0S17 
L 0287 

0 

.0384 

.0542 

,0839 . 

,1091 

.1329 

, 1565 

.1791 

.2012 

. 2246 

. 2472 

.2555 

.2048 

.0461 

.2070 

,4109 

.5308 

.6201 

.6837 

.7657 

.8402 

.9165 

.9914 

.9877 

0 

.0434 
,0847 
.1844 
' .1799 
,2266 
.2770 
.3339 
.3999 
44864 
.6065 
.7659 
' .9100 
.9563 
.7014 
.3249 
.0214 
.1992 
.3670 
.5063 
.64*4 
.7775 
.9012 
.9467 

1.8891 
. 5432 
.5435 
. .6859 
.5948 
.6008 
,-6086 
.6148 
.6210 
« 6314 
.6405 
, 6197 
. 4091 
.1212 
,8886 
.7151 
.8782 
.9800 
1.0262 
1,0529 
L 1018 
1.1288 
1,1437 

hom 

1.3891 
.5682 
.5740 
,6365 
.6656 
.6944 
.7297 
.7696 
.8190 
L 8932 
.9998 
1, 1271 
1.1743 
L0314 
.5399 
.0208 
.3259 
.5531 
.0991 
.8040 
.9100 
.9898 
1,0635 
1*0223 

1:8891 
,6732 
.6045 
.0870 
.7363 
.7888 
.8508 
.9244 
1. 0182 
L 1650 
1.3590 
L 6374 
1.8795 
1.9417 
L 4484 
.7567 
.2204 
.1262 
.3730 
.6651 
.7182 
.8509 
.9632 
.9013 

1.8891 
.5883 
.6850 
,7375. 
.8071 
,8817 
.9719 
L 0792 
1.2169 
L4168 
L 7183 
2,1477 
2.5848 
2.8620 
2.3568 
1.4925 
.7786 
.3007 
.0468 
.8062 
.5264 
.7119 
.8730 
.9403 

2.5600 
.9985 
1.0084 
1.0900 
1,1190 
1-1483 
1.1889 
L2378 
1.3038 
1.4090 
1,5697 
L. 7717 
1,8687 
1,7055 
1.0150 
.2607 
.2454 
.5622 
.7571 
.8852 
,999$ 
1.0739 
1.1201 
L 0580 

2.5500 
1.0185 
L 0340 
L 1407 
1.1897 
1.3419 
1:8100 
1J3926 
1.5024 
1<6708 
1.9289 
2.2819 
2.5789 
2.6157 
1.9234 
.9966 
.3069 
.1353 
.4310 
.6363 
.8080 
,9?49 
LQ299 
1- 0170 

2.5500 
1,0286 
1.0645 
1*1910 
1.2605 
1,3355 
1.4311 
1.5473 
1. 7010 
Li 9326 
2.2883 
2.7923 
3.2792 
3.5261 
2.8319 
1.7325 
.8592 
.2916 
.1049 
.8874 
,6162 
.7960 
,0300 
,9700 

2.8500 
1.0430 
1.0960 
12410 
1.8313 
14292 
16623 
1.7021 
18997 
21944 
20474 
3.8030 
3. 9844 
4.1864 
3,7403 
24083 
14114 
.7185 
,2213 
.1385 
.4343 
.0670 
.8494 
.9850 


Obtained by graphical Integration, 
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'ABLE III — ORDINATES AND VELOCITY DIS- 
TRIBUTIONS FOR NONSYMMETRICAL DUCT 
INLET WITH a=0.132 AND m= 1.5 

[t~ - 0.0341; * H »109°; * a =289.31°; * s -0,51°; * c -353.12°; (1=2,8687] 




at 

(ordinate) 


Velocity, o 

(C.-0.5; B= -0.38907) 


ci“0 

0,3 

0.6 

A-0.01818 

0.38901 

0.76285 

duot-lnlet seotlon 

L0000 

1.0000 

1.0000 

1,0751 

1. 1162 

1.1671 

1. 2487 

L3510 

1,4533 

1,2072 

1,3709 

L 6346 

1,1023 

1,3316 

1.6608 

,0626 

1.2689 

1. 5763 

.7630 

1. 1475 

1. 5712 

,4219 

.8743 

1,3207 

.0418 

.4934 

.9461 

,2706 

.1409 

.6614 

.5018 

.1377 

.2264 

. .644Q 

.3338 

.0236 

,7169 

.4693 

.2018 

.7479 

.5367 

.3236 

.7296 

.6600 

.3022 

,0994 

.6664 

.4316 

,6740 

.5663 

.4686 

.6661 

.6687 

.4813 

.6406 

.5700 

.4996 

.6262 

.6693 

.5134 

.6128 

.5697 

.6267 

,6986 

.5073 

.5361 

.5615 

.6324 

.6133 

,6469 

.5374 

.6279 

1,3891 

1.3891 

1.3891 




lower duot-lnlet seotlon 



TABLE IV— ORDINATES AND VELOCITY DIS- 
TRIBUTIONS FOR NONSYMMETRICAL DUCT 
INLET WITH «-0.2L5 AND m=2.0 

[r= -0.0374; * H ~114.98°; * 0 =289.I2°; *,,=6.19°; * 0 =383.12°; 6=2,8625] 


Velocity, v 

(r„~0.5; B- -0.3890?) 





Cl— 0 

A =0,03088 


Upper duot-lnlet section 


TABLE V— ORDINATES OF NON- 
SYMMETRICAL DUCT INLET 
WITH s=0. 280 AND m= 3. 0 
[b-1.6; t-0.0810; A— 3.00131 


Upper duot Inlet section 



0.1206 

1.0000 

.7002 

1,0666 

.9997 

1. 3818 

,7522 

1,2624 

.0176 

.9326 

,6288 

.6066 

.4788 

,1917 

.4427 

, 1627 

.4107 

,3887 

.3964 

.6611 

.3807 

.6669 

.3687 

.7284 

.3607 

,7660 

.3660 

.7797 

.3540 

,7663 

,3634 

.7217 

,3627 

,6929 

.8640 

.6717 

.3664 

.6631 

.3680 

.6360 

.3607 

.6203 

,3027 

.6039 

,3674 

.6660 

.8707 

.6436 

.1206 

1.3891 



Lower duot Inlet section 


Lower duot-lnlet section 



180.000 

CO 

191.230 

0.9017 

202.342 

.7604 

213.227 

.6618 

223.793 

.6799 

233.909 

.6222 

243.707 

.4817 

262.982 

.4473 

281.787 

.4232 

270. 130 

,.4024 

278.031 

,3895 

286,617 

,8802 

292,620 

.8799 

299.378 

,3826 

306.816 

,3940 

311,979 

.4119 

317.897 

.4464 

328,002 

.4926 

329.128 

.8787 

334,497 

,7170 

339.744 

.0830 

344.892 

1.8028 

849.908 

2,9347 

864. 906 

10, 1469 




VELOCITY DISTRIBUTIONS ON TWO-DIMENSIONAL WING-DUCT INLETS BY CONFORMAL MAPPING 




























































